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Abstract 

bo" 

We present an evaluation of the atmospheric tau neutrino flux in the 
energy range between 10 2 and fO 6 GeV. The main source of tau neutrinos 
is from charmed particle production and decay. The v T N — > tX event rate 
^ | \ for a detector with a water equivalent volume of 1 km 3 is on the order of 

\q \ 60-100 events per year for E T > 100 GeV, reducing to 18 events above 1 TeV. 

Event rates for atmospheric muon neutrino oscillations to tau neutrinos are 
■ also evaluated. 

OO 

\o 

<N ! I. INTRODUCTION 

t— i ; 

Recent measurements of the atmospheric neutrino flux by the Super-Kamiokande Collab- 



oration |0]|| show a deficit of muon neutrinos in comparison to theoretical predictions, while 
the measured electron neutrino flux is consistent with theory assuming that all neutrino 
masses vanish. Earlier, lower statistics experiments already showed inconsistencies with 



\ theoretical atmospheric flux predictions pi. On the basis of event rates and the zenith angle 
dependence of the muon neutrino deficit, the Super-Kamiokande Collaboration has shown 
that their results could be explained by neutrino oscillations between and u T 0. Oscil- 
lations imply at least one non-zero neutrino mass. Definitive evidence of massive neutrinos 
requires modifying the standard model of electroweak interactions. 

The Super-Kamiokande Collaboration measures neutrino fluxes from observations of elec- 
trons and muons in neutrino nucleon interactions: fi + N — > l+X. In view of the importance 
of the question of whether or not neutrinos have mass, one would like to see not just 
disappearance, but also the v T appearance coming from — > v T oscillations. Oscillation 
sources of z/ r 's include oscillations on the terrestrial scale from atmospheric u^'s as well as 
oscillations over large astronomical distances of v^s produced in, for example, active galactic 
nuclei [|J. A background to the flux of neutrinos from — > v T are tau neutrinos produced 
directly in the atmosphere. 

Tau neutrinos are produced in the atmosphere by cosmic ray collisions with nuclei in 
the atmosphere, which produce charm quark pairs. A fraction of the time, the emerging 
hadrons are -D s 's, which have a leptonic decay channel D s — > tv t with a branching ratio of 
a few percent. The subsequent r decays also contribute to the atmospheric v T flux. Heavier 
mesons contribute to the flux of tau neutrinos, but as we show below, they are negligible 
compared to the D s contribution. 
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In this letter, we outline the procedure to calculate the atmospheric tau neutrino flux. 
The details of the method as applied to atmospheric electron neutrino, muon neutrino and 
muon fluxes from charm decay appear in Refs. || and ||. We present our flux results for 
the neutrino energy range of 10 2 — 10 6 GeV, followed by the resulting v T + iV —> r + X event 
rates. For tau energies above 100 GeV, the rate is on the order of 60 — 100 events per year per 
km 3 water equivalent volume. With a 1 TeV threshold, there are on the order of 20 events. 
We also evaluate the expected event rate for the tau neutrino flux coming from — > v T 
oscillations based on a range of parameters consistent with the Super-Kamiokande results 
0. For tau energies above a few hundred GeV, the atmospheric tau neutrino background 
flux from D s — > v t t dominates the tau neutrino flux from atmospheric muon neutrino 
oscillations. 



II. TAU NEUTRINO FLUX CALCULATION 

The main source of atmospheric tau neutrinos is the leptonic decay of the D s : D s — > tu t , 
followed by r —> v T X. For relativistic particles, a semianalytic, unidimensional approximate 
solution to cascade equations describing proton, meson and lepton fluxes is a reliable ap- 
proximation [RBIfl. The solutions rely on factorizing source terms in the cascade equations 
into factors which are weakly dependent on energy times the incident cosmic ray flux, here 
approximated by a proton flux. The source term for pAir^ D s , for a D s of energy E and 
column depth X as measured from the top of the atmosphere is 



^ P (E,X) r dR p (E p ,O) \{E) dn p ^ D s 
\(E) Je p (f) p (E, 0) \{E P ) dE 



S(p-+D a )~ ™,'~ J l_ dE *^2 C)Z\ Z s E P ) (2.1) 



= ^ p JeY ZpdAE) - 

Here <j) p (E,X) is the flux of cosmic ray protons at column depth X. At the top of the 
atmosphere (X = 0), following Ref. 0, we set 

<f) p {E, 0) = 1.7 {E/GeV)- 2 - 7 cm-^^sr^GeV" 1 , (2.2) 

valid for energies values lower than 5 ■ 10 6 GeV. In Eq. (2.1), X p is the proton interaction 
length and dn/dE is the cross section normalized energy distribution of the D s emerging 
from the proton-Air collision. The quantity Z p D a (E) is called a Z- moment. Generically, Z- 
moments describe sources of particles of energy E, whether by production, decay or energy 
loss through scattering. A complete discussion of the Z-moment method of solution appears 
in Refs. and || . Its recent application to atmospheric muon, muon neutrino and electron 
neutrino fluxes from charm decays is found Refs. and ||. 

Solutions to the cascade equations in terms of Z-moments have two separate forms, 
depending on whether the decay lengths are short compared to the height of production 
("low energy") or long ("high energy"). For 2/ r 's from D s (and r) decays, we confine our 
attention to the neutrino energy range 10 2 — 10 6 GeV. These neutrino energies are well 
below the critical energy of ~ 10 8 GeV, above which decay lengths of the relativistic -D s 's 
and t's are longer than the vertical distance to height of production. Tau neutrinos are 
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called "prompt" in the low energy regime. The approximate solution for the v T + v T flux, at 
the surface of the Earth, is 



Kk e ) = — x _ z — • ( 2 - 3 ) 

The prompt tau neutrino flux is isotropic. Z pp (E) accounts for the proton energy loss in 
proton-air collisions. For Z pp (E), we use the results obtained by Thunman et al. in their 
recent evaluation || using the Monte Carlo PYTHIA |J. A similar, energy independent 
value was used in Ref. |§. 

For Z p d s (E), there are several approaches. Here we show the results from next-to- leading 
order (NLO) perturbative QCD production of charmed quark pairs, scaled by a factor of 
0.13 to account for the fraction of c — > D s (K|. Details of the NLO calculation in the 
context of the prompt muon, muon neutrino and electron neutrino fluxes appear in Ref. || . 
A second evaluation relies on Thunman et al.'s Z pD o(E), rescaled by the ratio of D s to 
D° production taken to be 0.25. 

The D s — > v T decay Z-moments have several contributions. The most straightforward is 
the direct D s — > v T in the two body decay, where |7]|| 



z {2 body) (E) = i 1 -** dx Z pDa {E/x)) a pA (E/x) <j> p {E/x, 0) B 
DsUt 1 ' Jo x Z pDs (E) cr pA (E) <f> p (E,0) 1-R D . 

in terms of x — E/Ed s , Rd s = rnl/rn 2 D s and B = 0.043, the branching ratio for D s — > tv t 
[ TT]] . We use the inelastic proton-air cross section 



a pA (E) = (290 - 8.71n(£/GeV) + 1.141n 2 (£/GeV)) mb (2.5) 



parameterized in Ref. JT3 

For neutrinos produced in the chain decay D s —> r — > v r the decay Z-moment is given 

by 

f 1 dy r 1 dx Z pDs (E/(xy))a pA {E/(xy))4> p (E/(xy),0) 
DsVt 1 ' Jo y Jr Ds x Z pDs (E) a pA (E) 4> P ( E ,°) ' 
B dn T ^,, 



1 - R Ds dy 



where B, Rd s and a pA are the same as in the previous case while now x = E t /Ed s and 
y = E/E T . The y distribution that appears in Eq. (2.6) can be generically written in the 
following way 

ihh -"- -- B T (g (y) - P T (x) gi (y)) (2.7) 



dy 

where B T is the branching ratio for r — > u T X and P T (x) is the tau polarization. We have 
evaluated the functions go(y) and gi(y) for r — > u T p, t — > u t tt, t — > u T ai [0 and r — > v r lv(_. 
They are collected in Table I, as are the branching fractions that we use. The purely leptonic 
decay also appears in Refs. In terms of the energy of the parent D s (Ed s ), the tau 

polarization can be written as 
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p = 2Rds E Ds (1 + R Da ) ^ ^ 

1 — RD a E t (1 — Rd s ) 



The total D s —> v T Z-moment is 



rv rvVl body) ryichain) / n n \ 

Ds^t = Z D 3 u T + T) s v T ) I 2 - 9 ) 



7 (2 body) . y(chain) 

^D a v T — 

7 (chain) 



where Z K D v includes the sum over all the tau decay modes in Table I. 



III. TAU NEUTRINO FLUX RESULTS 

Using NLO perturbative QCD with m c = 1.3 GeV, and factorization scale M equal to 
twice the renormalization scale \i = m c with the CTEQ3 parton distribution functions ||14|| , 
we obtain Z p d s between 10 -5 and 10 -4 . The Thunman et al. PYTHIA results (TIG) 
give instead Z v d s between 3.5 ■ 1CT 5 and 6.5 • 10~ 5 . The corresponding decay moments Zd sVt 
are 10 -3 — 5.5 • 10~ 3 and 0.8 • 10" 3 — 4.8 • 1CT 3 , respectively, over the energy range 10 2 — 10 6 
GeV. 

The results for the tau neutrino fluxes scaled by E 3 are shown in Fig. 1 where the 
solid line represents the flux obtained using NLO perturbative QCD, and the dashed line 
represents the flux obtained using the TIG D s production Z-moment. In principle, one 
should also take into account the neutrinos produced by the semileptonic decay of the b 
quark (b — > crv T ) followed by the decay of the tau, but this contribution is less than a 
few percent of charmed meson source. The dotted line represents the flux from b quark 
production and decay in the limit of zero polarization for the tau. 

The fact that the flux calculated using NLO perturbative QCD differs from the TIG 
PYTHIA based calculation || is apparently due to the inclusion of fragmentation effects in 
the PYTHIA Monte Carlo. A detailed comparison of the two approaches appears in Ref. 
0]. The tau neutrino flux evaluated using NLO perturbative QCD can be parametrized as 

log lo (E 3 i , T ( J E;)/(GeV 2 /cm 2 s sr)) = -A + Bx - Cx 2 - Dx 3 (3.1) 

where x = log 10 (£/GeV), A = 7.08, B = 0.765, C = 0.00346 and D = 0.00349. The TIG 
rescaled result calculated here has A = 6.69, B = 1.05, C = 0.150 and D = -0.00820. 



IV. DISCUSSION 

Large underground experiments may detect z/ r -nucleon charged current interactions v T + 



N —> r + X. In the absence of v u — > v T oscillations, event rates [15] for E T > 100 GeV and 
E T > 1 TeV are shown in Table II. We have summed over neutrino and antineutrinos and 
assumed that there is no attenuation of the v T flux due to passage through the Earth at 



these energies ||15|| . Detection of atmospheric v T 's in the absence of oscillations requires large 
detector volumes, on the order of 1 km 3 water equivalent volume for ~ 60 — 100 events per 
year with a tau energy threshold of 100 GeV. Such large detector volumes for the detection 
of u u charged current interactions are being considered [TO. 

In the case of neutrino oscillations, atmospheric muon neutrinos may oscillate as they 
travel to the detector. Even with a small fraction of v a — > v T conversions, the copious 
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production of u^s in the atmosphere make this a potentially large source of iv's. Again, 
assuming no attenuation due to neutrino interactions in the Earth, and approximating the 
height of production in the atmosphere by a constant 20 km ||17||, we have evaluated the 



probability of oscillation in terms of neutrino energy and zenith angle. Using the Bartol 
muon neutrino flux |18[ from tt and K decays extrapolated to E = 10 6 GeV, adding the 



NLO prompt flux || and integrating over zenith angle, Table III shows the event rates 
for two limiting values of Am 2 , the difference between the squares of the tau neutrino and 
muon neutrino masses, quoted by the Super-Kamiokande Collaboration 0. We have set 
the oscillation mixing angle sin 2 (26 l ) = 1. For the larger mass difference, — > u T events 
overwhelm the prompt v T events for E > 100 GeV, while for the lower mass difference 
the rates are comparable. Tau neutrino oscillations to in this parameter range have a 
negligible depletion of the prompt v T — > r event rate. 

For E T > 1 TeV, because prompt fluxes decrease less rapidly with energy than the muon 
neutrino flux coupled with oscillations, the prompt z/ r 's dominate the event rate. The cross 
over to prompt v T domination of the flux occurs at a few hundred GeV. Another distinction 
between prompt and oscillation tau neutrinos is the angular distribution. The prompt source 
is isotropic, while the oscillation to z/ T 's will be strongly angular dependent, with the most 
zz/s traveling upward because of the long path length. 

We have concentrated here on atmospheric sources of prompt z/ r 's and atmospheric — > 
v T . There are many flux models of astrophysical sources of z/ M 's ||]. If neutrinos have mass, 
because of the large astronomical distance scales, on the order of half of the z/^'s would 
oscillate to z/ r 's. The detectability of v T — > r conversions in the PeV energy range has been 
explored in the literature |19|]. It is clear from our calculation that the atmospheric tau 



neutrino flux is negligible at PeV energies. 
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TABLES 



TABLE I. Functions go and g\ in the tau neutrino energy distribution from r decays, in terms 
of y = E/E T and T{ = mf/m 2 . and relative branching ratios. 

Process B T go g\ 

T^VrltVp 0.18 |_ 32/ 2 + |y 3 l_ 3y 2 + 8 y 3 



u T p 0.26 



2£-l+rp \ / l-2r p 
l-r p I (l_r p )* J I l+2r p 



n iq 1 I 2y-l+r ai \ / l-2r ai 



TABLE II. Charged current event rate per year per km 3 water equivalent volume from the 
prompt u T + v T flux. 

Threshold NLO QCD TIG 

100 GeV 58 98~ 

1 TeV 18 18 



TABLE III. Charged current event rate per year per km 3 water equivalent volume from 
v ii + —> v T + v T oscillations, assuming sin 2 (2#) = 1. 



Threshold 


Am 2 = 5 • 10~ 4 eV 2 


Am 2 = 6 • 10~ 3 eV 2 


100 GeV 


71 


9100 


1 TeV 


0.036 


5.2 
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FIG. 1. Tau neutrino plus antineutrino fluxes from D s using NLO perturbative QCD (solid 
line) and the TIG rescaled Z- moment (dashed line). The tau neutrino plus antineutrino fluxes 
from b quark NLO perturbative QCD production and decay with the tau polarization equal to zero 
is shown by the dotted line. 
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